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ABSTRACT 

A mathematical model taking into 
account small (and constant) 
gravitational levels is developed for 
vaporization of an isolated liquid 
droplet suspended in a stagnant 
atmosphere. A goal of the present 
analysis is to see how small 
gravitational levels affect droplet 
gasification characteristics. Attention 
is focused upon determining the 
effects on gas-phase phenomena. The 
conservation equations are normalized 
and nondimensionalized, and a small 
parameter that accounts for the effects 
of gravity is identified. This parameter 
is the square of the inverse of a Froude 
number based on the gravitational 
acceleration, the droplet radius, and a 
characteristic gas-phase velocity at the 
droplet surface. Asymptotic analyses 
are developed in terms of this 
parameter. 

In the analyses, different spatial 
regions are identified. Near a droplet, 
gravitational effects are negligible in 
the first approximation, and the 
flowfieid is spherically symmetric to 
the leading order. Analysis shows, 
however, that outer zones exist where 
gravitational effects cannot be 
neglected; it is expected that a 
stagnation point will be present in an 
outer zone that is not present when 
gravity is totally absent. The leading- 
order and higher-order differential 
equations for each zone are derived 
and solved. The solutions allow the 
effects of gravity on vaporization rates 
and temperature, velocity and species 


fields to be determined. 

NOMENCLATURE 

D mass diffusion coefficient 
e internal energy 

<=Sh,Y, - P/p) 

1*1 

«* 

f normalized gravitational vector 
F r Froude number (= U, / -^gR) 

g gravitational acceleration 
L heat of vaporization 

Le Lewis number (= X g /(p g C p D g )) 

Ma Newtonian Mach number 
m evaporation rate 
P gas-phase pressure 
R instantaneous droplet radius 
Re Reynolds number (= U.R / v_ 

o' 

SR gas constant 

Sc Schmidt number (= v„ /D ) 

g/°° ' g/ oo/ 

u gas-phase velocity 
T gas-phase temperature 
U a radial velocity at the surface of 

the droplet (= — — ) 

Pg 8R' 

Vi diffusion velocity of species i 
Y| mass fraction of species i 

z contracted radial coordinate 

(= e^r) 

T] contracted radial coordinate 
(= e^r) 

m 

x stress tensor 
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e Richardson number (= gR/U s ) 

<J> viscous dissipation of 
energy 

y ratio of specific heats 

(=C p /C v ) 

v kinematic viscosity 
p. dynamic viscosity 
p gas-phase density 

Subscripts and Superscripts 

* dimensional quantity 
~ outer expansion variable (zone II) 

A outer expansion variable (zone III) 

g gas phase 

d liquid phase 

r radial direction 

s droplet surface 

0 leading order quantity 

1 higher order quantity 

OO far field state 


INTRODUCTION 

The last few decades have seen an 
increasing interest in theory and 
experiments on evaporation and 
combustion of a liquid droplet. This 
interest reflects the importance of this 
area to the understanding of the 
fundamentals of combustion science. It 
also reflects the fact that a significant 
portion of current practical 
applications use liquid fuels as a source 
of energy due to convenience in 
transportation and storage. In most of 
these applications spray combustion is 
used, where liquid fuel would exist as 
droplets inside a combustion chamber. 

A detailed study of spray combustion is 
difficult due to the complexity of 
interactions between droplets. One 
logical way to gain an understanding 
of the general behavior of spray 
phenomena is to study droplet 
vaporization and/or combustion under 
well-controlled conditions. When 
experiment and theory agree for 


simplified situations, predictions may 
be made with increased confidence for 
more complex cases where accurate 
experimental data may not exist. For 
this reason, many studies have focused 
upon isolated droplets undergoing 
spherically symmetrical vaporization 
with and without combustion. 

Reducing the acceleration of gravity to 
negligible levels is one of the most 
critical conditions for isolated droplet 
experiments that strive for spherical 
symmetry. Based on this condition, the 
spherically symmetric vaporization 
and combustion assumption can be 
applied provided that there is no forced 
convection or any interaction between 
the droplet and its environment other 
than heat and mass transfer. However, 
reducing the gravitational 
acceleration is not a simple task and a 
nonzero effective gravity level is 
invariably present in these type of 
experiments. The effects of small 
gravitational levels on droplet 
phenomena are largely unknown. 
Therefore, studies to model the effects 
of gravity on the vaporization and 
combustion of isolated droplets are 
warranted. These studies can increase 
understanding of the importance of 
small buoyancy levels and in turn will 
help in the interpretation of results 
from microgravity experiments. 

In this study, we will develop a 
mathematical model using asymptotic 
expansions for an isolated liquid 
droplet suspended in a stagnant 
oxidizing atmosphere which undergoes 
diffusion dominated vaporization, 
taking into account gravity effects. In 
many recent theoretical droplet 
models, buoyancy effects have been 
neglected. In contrast, gravity levels 
are never zero in experiments. For 
example, the effective gravity level 
can be reduced to 10*® m/s 2 (at best) in 
the space shuttle (though accelerations 

of about 10* 3 m/s 2 are more typical), to 
10* 5 m/s 2 in drop-towers, and to 10* 3 
m/s 2 in a drop-tube apparatus 1 (Wang. 
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D. and Shaw, B. D., 1991). Depending on 
the droplet diameter these gravity 
levels can be important and might 
significantly affect outcomes of the 
experiments. There are quite a few 
published articles that cover various 
aspect of droplet evaporation and 
combustion, but there is a no published 
source that looked at the buoyancy 
effects within this frame work. In this 
study, we hope to shed light on this 
factor and provide a basic 
understanding of its importance. 

THE GOVERNING EQUATIONS 

The dimensional governing equations 
can be written as follows : 

3p / dt + V . (p U ) = 0 ( continuity ) 

3Yi/8t* + u* . VYi 
= [^’(p Y i^i)]/ P* ( spedes ) 

p*3e*/3t* + p*u*. Ve* 

= V.( VVT* ) - P* ( V. U* ) (energy) 

p*9u* / 3t* + p*u*. Vu* 

= " + V .X + ^p* - pi, jg* (momentum) 

P = P 91T* ( equation of state ) 

where: 

hi* = h|*° + J a pA dT* 

T° 

Let VO = p^. g* 

P*=p*-<D 


These equations will be 
nondimensionalized by using the 
following variables : 



where U s * is the radial velocity at the 
surface of a droplet undergoing purely 
spherically symmetrical combustion in 
the absence of gravity. 

ASSUMPTIONS 

For a liquid fuel droplet vaporizing in an 
oxidizing atmosphere of an insoluble 
gaseous species, the above equations are 
simplified according to the following 
assumptions : 

Quasi-steady state vaporization, so the 
time-dependent terms can be neglected. 
Taking the molecular weights and the 
specific heats to be constants, and the 
product pD to be constant. Constant gas- 
phase transport properties. Lewis and 
Schmidt numbers are constants. Perfect 
gas. Single fuel species droplet. No spray 
effects, that is the liquid fuel droplet is 
isolated and suspended in a stagnant 
environment. Constant thermodynamic 
pressure. Constant and uniform droplet 
temperature. The Newtonian Mach 
number is small, M 2 « 1. Fick's law for 
mass diffusion ( V 4 = -DVlnYj) holds. 

Using these assumptions as well as the 
nondimensional variables allows the 
governing equations to be 
nondimensionalized as follows: 

V . <pu> = 0 (continuity) 
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ANALYSIS TECHNIQUE 


Re Sc pu . VYt 

= V 2 Yj (species) 

Re Sc pu . VT - Le V 2 T 

= (y - 1)M 2 (({) + U. Vp) (energy) 


pu . Vu + 
V.t 

“ Re ^ 
P= pT 


vp 

Ma 2 

e(p-l) f 

(state) 


(momentum) 


The dimensionless boundary conditions 
are : 

at r = 1 
U0 =0 
T = T S 


N 

XY, = 1 


1=1 


Y = 



y l 

y-1 



M 


J 


8Y, 



Re Sc u r 

5 


(Y 


is 


- 1 ) 


Re Sc u r L 
LeD 


and as r -» 
u r — > 0 
u 0 -* 0 
Y, -» Yu 
T 1 


oo 


p -> 1 
VP 
Ma 2 


e(p-l)f 


The technique of matched asymptotic 
expansions will be used as a tool to 
study the effects of gravity on the flow 
variables during the evaporation of a 
liquid fuel droplet. 

There have been a number of studies of 
droplet evaporation and combustion 
using this perturbation technique. For 
example, Fendell, Spankle and Dodson^ 
(1966) used a singular perturbation 
approach to consider droplet 
combustion with convective flow field 
that is entirely in the Stokes regime. A 
Reynolds number was used as their 
small parameter of expansion. Fendell, 
Coats and Smith also studied droplet 
vaporization in a slow, compressible 
and viscous flow^. Kassoy et al.* (1966), 
studied compressible low Reynolds 
number flow around a sphere. They 
derived inner and outer expansions for 
the flow variables in terms of the 
temperature difference, Reynolds 
number and the free stream Mach 
number. Acrivos and Taylor 5 derived 
an expansion in terms of a small 
temperature difference between a 
sphere and free stream and calculated 
the average Nusselt number of the 
incompressible flow over a sphere 
with a Stokes velocity profile. Gogos et 
al.6 studied evaporation of a fuel 
droplet with a strong evaporation- 
induced radial velocity while 
undergoing slow translation. They 
introduced the translation as a 
perturbation to an otherwise 
stationary droplet. Gogos et al. have 
also considered combustion of droplets 
undergoing slow, translation^. 
Mahoney®, Fendell^, and Sato 1 ® have 
presented asymptotic analyses of 
natural convection heat transfer from 
rigid spheres assuming that the 
Grashof number is small, while Hieber 
and Gebhart 1 1 studied mixed 
convection about rigid spheres in the 
limit of small gravity effects. 
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In the model presented here, a quasi- 
steady state assumption is used because 
of the fact that the ratio between the 
densities of the gas-phase and of the 
droplet is small 0(10" 3 ). Matched 
asymptotic expansions are used based 
on a small parameter identified as the 
ratio of the inertia force to the gravity 
force. This parameter is the inverse of 

the Froude number (U\/yg*R*) 

squared, where U% is the droplet 
surface ‘velocity, R is the instantaneous 
droplet radius, and g* is the gravity 

acceleration. This parameter is the 
variable e that appears in the 
dimensionless momentum equation. 

PHYSICAL CONSIDERATIONS 

An order-of-magnitude analysis of the 
governing equations indicated the 
relative importance of each term in 
the equations. It was found that the 
order of magnitude of the non- 
dimensional gravity term in the 
momentum equation was a useful guide 
to distinguish the different governing 
physical characteristics. This term was 
equal to the reciprocal of the Froude 
number squared and is denoted as 

e = g*R/U* 2 , where £«1. For typical 
situations, a 1 mm droplet has values of 
£ ranges from 10" 4 in drop-tube 

apparatus to KT 6 in space shuttle 
experiments. In order to study the flow 
field around a droplet, the governing 
equations have to be solved. These 
equations in their original form are 
very complex. Simplifying assumptions 
have to be made in order for one to 
approach these equations analytically. 
For instance; the thermodynamic 
pressure was considered to be constant 
since characteristic Mach numbers 
(M) for the flow are significantly less 
than unity. The pressure gradient in 
the momentum equation, on the other 
hand, can not be neglected since it is 
divided by (M 2 ) which is very small. 


The non-dimensional equation of state 
is then simply pT = 1. In the energy 
equation the second term of the right 
hand side can be neglected. This 
dissipation term is of 0(M 2 ), where M 
« 1. Taking the Lewis number as 

constant will allow us to have 
pD ^constant. 

A physical approach was used to deduce 
that this problem is singular for £ -4 0, 
as well as the expansions for inner and 
outer regions. In an inner region near 
the droplet (termed zone (I)), the first 
term in the expansion is the 
spherically symmetric solution. This 
solution will prevail in the absence of 
any significant gravity effects and the 
flowfleld will only vary in the radial 
direction. For the next term in the 
expansion, a higher-order correction 
of 0(E) will be applicable to the 
velocity field, which is a result of the 
body force term. From the continuity 
equation it can be seen that the order 
of magnitude change in density is as 
large as the change in velocity. From 
the equation of state, the order of 
magnitude changes in temperature 
will be as large as the changes in 
density. Therefore, it was concluded 
that the 0(£) correction is applicable 
to all of the flow variables (velocity, 
density, mass fraction, temperature 
and pressure). To locate zone (II) i.e., 
where gravity is important, we had to 
determine where the spherically- 
symmetrical solution breaks down. 
From the leading-order continuity 
equation, the velocity field was 
obtained. This velocity field was then 
used in the radial momentum equation 
to estimate where the largest term of 
the equation becomes of the same order 
as the gravity term, i.e. 0(£/r). It was 

found that the largest term of the 
inertial, pressure and viscous terms is 
of 0(l/r 5 ). Thus, it was concluded that 
the negligible buoyancy assumption is 

no longer valid when (1/r 4 ) « £, and 
that gravity should be accounted for 
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when r - 1/E 1 / 4 . The rescaled radial 

coordinate z = £^ 4 r was introduced and 
expansions for all flow variables were 
obtained. 


SOLUTION METHOD 

The method of matched asymptotic 
expansions is used in this study to solve 
the governing equations that describe 
the flow field around the droplet. 
Inner .and outer expansions are 
derived for each of the flow variables, 
(velocity, density, mass fraction, 
temperature and pressure), in terms of 
a small parameter (£). The solutions 
are typically developed up to the 
second terms of the expansions. The 
basic simplification of the governing 
equations comes from the assumption 
that £«1 where as stated above, the 
spherically symmetric assumption is 

only applicable when (l/r 4 )» 6, 

since when (1 / r*) =* £ the effect of 
buoyancy would be felt in the flow 
field and this will invalidate this 
assumption. 

In zone ( I ), the gravity term suggests 
an 0(E) correction to the velocity field. 
This correction is also applicable to the 
other dependent variables. 

The inner expansions valid near the 
droplet are as follows : 


u(r,0 ; e) = u 0 + zu x +. 

p(r,0 ; e) = p 0 + epi + 

Y(r,0 ; e) = Y 0 + eY 1 + . . . . 

T(r,0 ; e) = T 0 + eT, + 

P(r,0 ; e) = P 0 +■ ePi + 

where : 

Ui =[u r i O (r) + u rU (r)cos0]e r 
+ [u eil (r)sin0]e e 

Pi =Pio( r ) + Pn( r ) cos0 
Y, = Y lo (r) + Y n (r)cos0 
T| = T lo (r) + T n (r)cos0 
Pi =P lo (r)+P n cos0 

Substituting the above expansions into 
the governing equations produces the 
following equations: 

Order c®: 


V . (p 0 u 0 ) = 0 

Re Sc p 0 u 0 • VY 0 = V 2 Y 0 


Re Sc s_ 

— — p 0 u 0 . VT 0 = 
Le 

n„ . Vii. = --^Eo 
Ma 2 


p 0 u 0 . Vu 0 


V 2 T 0 

V .to 
+ Re 


P0T0 ~ 1 


Boundary conditions : 
atr = 1 



YL 

Y-l 




s 


TT 

b s 


J 


3Y 0 | = ReScu ro 

“drl, D 
aT 0 | = RfiScL u 
1TL LeD 
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The analytical solution to the above 
spherically symmetric ordinary 
differential equations is : 



Yq = 1 + (Y 0o - 1) expl-ReSc- 


-Ma 2 !'*— ° 2 T 2m o 2 dT 0 

dr Ma f r 5 T ° —IF 

+ 8 f m o d 2 T 0 m 0 dT 0 1 

3Re|_ r 2 dr 2 r 3 dr 
where m 0 is the integration constant 


Order e 1 : 

V . (p 0 u, + p t u 0 ) = 0 

Re Sc ((poit! + Pl u 0 ) . VY 0 
+ p 0 u 0 -VY! ) = V 2 Y! 

((PoGt + PtUo) • VT 0 
+ PoQo-VT 1 ) = V 2 T 1 


P 0 u 0 • Vu 1 + (p 0 u 1 + pju,,) . Vu 0 


VPt V.Ti 
Ma 2 Re 


+ (Po-i)? 


pi - - Tj 


The above equations will actually consist of 
two sets of equations: 


[ r2 (PO u rlO + PlO u rC>)] = 0 (continuity) 

m 0 du r iQ , m 10 dUrO ( 1 dP 1Q 

r 2 dr r 2 dr Ma 2 dr 

= + (momentum) 


Re[3 dr 2 3dr^ r 


Re Sc f m 0 dT 10 m 10 dTp ^ 

Le V r 2 dr r 2 dr J 

" 7 s ( r2 ^ a ) =0 <ene, ®’ , 

ReScf BL jlm + 5Uo 

V r 2 dr r 2 dr J 

1 df 2 dY 10 > \ 

~7to[ T "drj =0 <species) 


PlO = - Po T io 
atr = 1 
Y ilO = 0 


(state) 


Tin = 0 



Re Sc 


[ u rlo(l“ Y o) + u rO Y lo] 


LU . 


Re Sc L 
LeD 1 



Applying the method of the integrating 
factor to the energy and species equations, 
we get: 


at r — 1 , Uj-q — 1 

in . r(Pod) u rlO^ + Plo( 1 )) 

■" *T 1+ T Sc5 


Tio(r)= _^ (T ,_ L)+ ^ 

^m 10 f Re Sc f 


Re Sc, 


Le 
+C n expl 


jexp 

{ ' ReScmp A 
l Le rj 


Le 


m 0 



where C| and C ^ are constants 


V (r>\ — m !0 I ^*2 

10 m 0 ReScm 0 
-Re Sc -Dexp^-ReSc-S-j 

+C 22 exp^- Re Sc 

where C2 and C22 are constants 


Evaporation Rate Modification: 

Let |i = cos0 

Then the evaporation rate can be defined 
by: m = 27cr 2 J pu r d|i , 

substituting the inner expansions into 
the above equation, we get : 

m = 2rcr 2 [2pou r o + 2e(ppu r ip + Pio u ro)+— ] 

comparing the above equation to the 
spherically symmetric rate of 

evaporation, m = 4jt(r 2 pou r o), we get : 


From continuity in the first order 
problem we have : 

= p 0 (r)u rl 0 (r) + Pi 0 (r)u r 0 (r) 
r 

=> ntjo = (po(l) u rlO^) + P10(D) 

Where mjo can be found from the 
solution of the energy or the species 
equations for r = L 

By substituting the second part of 
the higher order expansion we get 
coupled boundary value problems 
that can be solved numerically. 


drirll 

dr 


' = _f2 + _T_dp a '! 
V r Po dr ) 


. m Q dT n 
* 


UrU ' 2 dr 


+BLlU.dP!>_2HS 1 l =0 (continuity) 

PO r 2 dr 


d 2 uen _ fm 0 2^ duel 1 1 du rll 

dr 2 i, r 2 r J dr r dr 
Pll , 4^ dUfii Urii+uen Y 


<■ 


Ma 2 3 1 dr 
m 0 4 'l 


j. 


u rll 


+ pp (0 - momentum) 


m_ = u c r (Po(r)Urlo(r) + Pio(r)u r0 (r)) 

mp L po(r)u r0 (r) 
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a ll 


F 11 | 4 T du r ll u r11 +u eil 
Ma 2 3 v dr r 


^ a ii _ f m o 4^ du r n 2 du e n 
dr v r 2 rj dr r dr 
( du. 

{*>-£ 


r 2 y 


u rll + — 2' u 011 


( — mg | du r n 

A 11 — 2~J— Pq (r - momentum) 

d2T ll _ ( 2 ReScm 0 _ V 
dr 2 "V "Lr r 2 Po J Tl1 
, ( Re Sc m 0 2 ^ dT^ 




Le 


J dr 


r 

dT 0 


, ReSc 

+ TT" rl,Po dr 
d 2 Y 


(energy) 






Re Sc 




2W 


n 


x) dr 
dY 0 


+ReScU rll p 0 — — — (spedes) 
dr 


Numerical solutions of these equations 
will be reported at a later time. 

In zone ( II ), the expansions for the 
flow variables will be of the following 
form: 

u(z,0;e) = £P n u n (z,e;e) 

n=0 

Pfe^e) = la n p n (z,0;e) 
n=0 

Y(z,0;e) = £8 n Y n (z,0;e) 

n=0 

T(Zr0;e) = SXnt n (z,0;e) 

n=0 

P( Z/ 0;e) = STi n P n ( 2/ 0;e) 
n=0 

where the coefficients P n , a n/ 8 n , 


% n , and T] n are, in general, functions 
of £. To determine the first coefficient 
of the above expansion, the zero th - 
order solution of the inner expansion 
was substituted into the radial 
component of the momentum equation. 
It was found that the viscous, pressure 
and convective-acceleration terms 
were 0(l/r5), while buoyancy effects 
were O(^r). Therefore, the buoyancy 
effects are negligible as long as 
(l/r 4 )» e; when (1/r 4 ) = e the 
above assumption is no longer valid. 
An outer variable (z) is introduced as a 

contracted coordinate, where z = 6^r. 
In the leading-order momentum 
equation, the density appears as a 
constant. In a constant-density flow 
field, the buoyancy will have no 
effects on the velocity field and 
therefore, the leading-order flowfield 
is spherically symmetrical. 

To determine the appropriate higher 
order corrections for the dependent 
variables in zone (II), the leading- 
order composite solution was 
substituted into the governing 
equations. The terms were then 
expanded. The momentum equation was 
found to be satisfied to leading order, 

with an error of O(e^) . This suggested 
that the higher order correction terms 
for the flow variables in zone (II) are 

0(e^) smaller than the leading-order 
terms. Thus, the coefficients then are : 

Po = # ct 0 = 8 0 = xo = 1 , Tl 0 = e 

Pi = / a, = 8, = xi = , ti 1 = e 

And the outer expansions which are 
valid far away from the droplet arc : 


9 



z = e^r 
V = 

ti(z,0; e) = e%uo + +..•••• 

p (z,0; e) = 1 + e^po + e%pi + 

Y(z,0; e) = Yoo + e^Y 0 + e%Yi + 

T(z,0; e) = 1 + e^T o + e%T i + 

P(z,0; e) = ePo + e%Pi + 

Equations for the correction terms 
were derived by substituting the outer 
expansions into the governing 
equations and then grouping terms of 
equal order of magnitude. 



Solutions to the energy and species 
equations are : 

To= i^-P.fcosfl) 

n*0 Z 

Y„ « £ Pn<“s8) 

n= 0 Z 

where P n is Legendre polynomial 
of order n. 

Matching the leading order solutions: 
Inner leading - order solution as r — > °<> 


Zone (II) leading-order equations: 
=0 

V 2 T 0 = 0 

f 2 Y 0 = 0 


u r ~^+ 0 (r- 3 ) 

r z 

T - 1 - ^^[1 - (T, - L)] + 0(r -2 ) 
Y ~ Y„ - (Y„ - 1 ) ReS ^ m ° + 0 (r - 2 ) 


v . T 0 = 0 
Po = -To 

as r — » oo 

u z o 0 

Ueo 0 
Tq — > 0 
Y o -*0 

In the leading-order equations, the 
solution is spherically symmetric and 
the effects of the buoyancy are not yet 
apparent. Only the radial velocity 
component is nonzero; this component 
is shown below: 


p _ 1+ ^Sci5a[ 1 _(T 8 -L)] + 0(r- 2 ) 

Le r L 

+CKr' 5 ) 

Zone n leading - order solutions as z -» 0 


U( iho^j 
u z ~ e 75 ! -y- 1 


T~l+e 


Y~ Y ~ + e ^[~r) + 
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Matching the two solutions gives : 
N n = S n = 0 for n £ 1 
m 0 =m 0 

XT Re Sc m 0 r , . 

N 0 = _o[i_( Ti _ L) ] 

V 0 = -ReSc m 0 (Y - - 1) 


The higher-order equations : 


V . (p 0 u 0 + Uj) = 0 

Re Sc _ =- = 2 - 

— u 0 .VT„=V 2 T, 

ReScu 0 .^Y 0 = ^ 2 Y x 

ao-Vu 0+ ^- = ^i + p n f 


Ma 


Re 


Pi - "Tj 


The continuity and the momentum 
equations can be solved using a form 
of a solution expressed in terms of 
the stream functions : 

u , - - 1 d V , m o N o 
z z 2 sin0 d0 z 3 
G _ 1 dy 

zsinO dz 

Substituting the above equations into the 
governing equations and taking the curl of the 
momentum equation yields the following: 

DV = sin 2 e^.J 

where 

D= d 2 + sine d f l d 
3z 2 z 2 dO^sinOdoJ 


as 


r — > oo 


u r3 0 
u 0 i -> 0 
Tj 0 

Y,-»0 

^■-6? 

Ma 2_Pof 

Solutions to the energy and species 
equations are : 


Ti = £ ^a-P„(cos9)-^*aNe 


n=0 Z 


U 2z 2 


*1 = X -%P„(a)se)-ReScm 0 ^2. 

n=0 z 2z z 

where P n is Legendre polynomial 
of order n. 


The general solution is found to be : 

V = f^ J -+ A 2 2 + A 3 z 2 + A,z 4 
+^ a -z 3 lsin 2 9 

a «i= 2 (^f- + '7- + A 3+A, z2 

+— zloosO+^^A 
4 J z * 

* ei = ~ ^ ~ 2 Aa " 4 A<z2 

-^zlsine 

8 J 

w ^ ere A|/ A2/A3 and are constants 


From the principle of minimum 
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singularity 1 2, it can be concluded that 
A4 should be equal to zero. The velocity 
profiles will then be as follows : 


u z i 


u ei 


-<■ 
-(*- 


M 3N 0 _~\ 


^. + A2 + No z ]cose + 51sNo 
j z 4 ' 


8 ) 


sin 9 


For the above solution, the first-order 
velocity is unbounded as z increases, 
while the leading-order velocity field 
decreases as 1/z 2 . Eventually, as z 
increases, the velocity correction term 
will become as large as the leading-order 
velocity term. When this happens, the 
solution in zone (II) breaks down, and 
more analysis is necessary. Another 
zone is necessary so that the solution can 
satisfy the outer boundary conditions. 


Zone (III): 

The new scaling factor in this zone is 
predicted from the failure of the 
solution of zone (II). the expansion 
equations for the this region are: 

X 

ti = e /3 r 

V= eXf 

u(t|, 8; e) = eX{i 0 + e uj + .... 

P (t|,8; e) = 1+ eXp 0 + eXpi+ 

Y(ti,0; e) = Yoo + eXy 0 + eXy 1+ 

- l/- V - 

T(ti, 0; e) =1 + £^3 Tq + . . . . 

P0n,9; e) = ePo + eXpj + 

The leading-order equations are as 
follows: 


V . (u 0 ) = 0 
^ 2 T 0 = 0 


A * A 

V 2 Y n 


= 0 


Ma 2 Re 
Po = "T 0 


as r 00 
“zo “♦0 

u 0O ->0 

T 0 ->0 

Y 0 -»0 

The solutions to the energy and species 
equations are : 

To= t^ncosB) 

n-OTl 

Y 0 = I^rPCcosG) 

n-O^l 

where P„ is Legendre polynomial of 
order n. 

Again, matching the leading - order 
solutions in zone I with the leading - 
order solutions in zone II, we get : 


D n = S n = 0 for n £ 1 

mg = m 0 = mg 

D 0 =N,=- Re ^ ro ° [l-(T,-L)] 

S 0 = V 0 =-ReScm 0 (Y.-l) 

The velocity profile can be obtained 
using the stream function approach : 

„ 1 dy 

1,0 T) 2 sin 9 d0 

. 1 d\pr 

U0 ° Tisin0 dri 

Substituting the above into the 
momentum equation and taking the 
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curl yields: 


DISCUSSION AND CONCLUSIONS 



where 


p _ 9 2 | sin9 d_f 1 <n 

dr] 2 ri 2 dslsinBdsJ 

The general solution is found to be : 

V = + B 2 TI + B 3 t| 2 + B 4 t| 4 + "H 3 jsin 2 0 


^0 = 2 ■^- + ^- + B 3 +B 4 ri 2 +^-Tijcos0 


ueo = 2B 3 -4B 4 Ti 2 -^&T|jsin0 

where Bj,B2,B 3 and B 4 are constants 


The nondimensionalization of the 
equations produced the result that the 
variable e characterizes the influence 

of gravity. When e is small, the 
flowfield is momentum-dominated near 
the surface of a droplet, with buoyancy 
subdominant. Since velocities decrease 
as fluid particles travel away from the 
droplet surface, buoyancy effects 
eventually become as important as 
momentum and viscous effects. 

For 6 — » 0, the analysis of zone (I) 
indicates that the flowfield near the 
droplet surface will be spherically- 
symmetrical to leading order. 

Correction terms are of order e in this 
zone for all variables. A result of this 
analysis is the expectation that droplet 
vaporization rates should increase 

linearly with e so long as 6 « 1. 
Eventually, however, as r increases, 
the solution in zone (I) breaks down, 
which leads to the analysis of zone (II). 


From the principle of the minimum 
singularity (see Van Dyke), it can be 
concluded that the constant B 4 should 

equal to zero. The velocity profiles arc as 
follows: 

,» ( Bi B 2 3Do ^ 

u 90 = 3 — —•‘H sine 

U 3 T] 8 J 

As was found in zone (II), the velocity 
field in zone (III) increases without 
bound as T) increases. Again, this 
implies that analysis of at least one 
other zone is required to complete the 
analysis of this problem. 


In zone (II), the leading-order 
solutions are spherically symmetrical, 
while the correction terms include 
nonzero velocities in the 0-direction 
that tend to turn the flow in the 
direction of the gravity vector. 
Eventually, as the coordinate z becomes 
large in one (II), the asymptotic 
correction term becomes as large in 
magnitude as the leading-order term. 
This is interpreted as an indication that 
a stagnation point will likely exist near 
the outer boundary of zone (II). To 
reinforce this idea of the existence of a 
stagnation point, inspection of the 
solutions obtained for zone (III) 
indicates that the velocities in this 
zone will in general be "downward," 
that is, they will be in the direction of 
the gravity vector; these downward 
velocities should produce a stagnation 
point along a symmetry line of the 
flowfield. Based upon the scalings 
derived in this paper, the 
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dimensionless radial location (r s ) of 
the stagnation point is expected to scale 
as r s * A qualitative sketch of the 

expected flowfield in the vicinity of a 
droplet is shown in Fig. (1). Only the 
streamlines on one side of the droplet 
are shown; the streamlines on the 

other side are symmetrical. 

As noted previously, further analysis is 
required to provide a description of the 
flows outside of zone (III); this 

research is presently underway and 

will be. reported at a later time. Based 
upon previous analyses®*^’ *•*., it is 
expected that a plume-like structure 
will eventually be formed below the 
droplet; boundary-layer theory can 

likely be used to analyze the plume 
structure. However, analysis of the 
transition from the flowfield near the 
droplet to the plume-like flow is 
expected to be challenging. A sketch of 
the expected flowfield far away from a 
droplet is shown in Fig. (2). Only the 
streamlines on one side of the droplet 
are shown; the streamlines on the 
other side are symmetrical. 
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